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We consider attractive irreducible conservative particle systems 
on Z, without necessarily nearest-neighbor jumps or explicit invariant 
measures. We prove that for such systems, the hydrodynamic limit 
under Euler time scaling exists and is given by the entropy solution 
to some scalar conservation law with Lipschitz-continuous flux. Our 
approach is a generalization of Bahadoran et al. [Stochastic Process. 
Appl. 99 (2002) 1-30], from which we relax the assumption that the 
process has explicit invariant measures. 

1. Introduction. In this paper we study the hydrodynamic behavior of 
a class of asymmetric particle systems of Z, which arise as a natural gen- 
eralization of the asymmetric exclusion process. For the latter a variety of 
results and approaches are available: the hydrodynamic limit is given by the 
entropy solutions to the scalar conservation law 



where u(-, •) is the density field, G — the macroscopic flux — is given by G(u) = 
7it(l — u), and 7 is the mean drift of a particle. Because there is a single 
conserved quantity (i.e., mass) for the particle system, and an ergodic equi- 
librium measure for each density value, (1) can be guessed through heuristic 
arguments if one takes for granted that the system is in local equilibrium. 
The macroscopic flux G is obtained by an equilibrium expectation of a mi- 
croscopic flux which can be written down explicitly from the dynamics. For 
the simple exclusion process, equilibrium measures are also explicit, and so 
is G. A rigorous proof of the hydrodynamic limit turns out to be a diffi- 
cult problem, mainly because of the nonexistence of strong solutions for (1) 
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and the nonuniqueness of weak solutions. Since the conservation law is not 
sufficient to pick a single solution, the so-called entropy weak solution must 
be characterized by additional properties; one must then look for related 
properties of the particle system to establish its convergence to the entropy 
solution. 

The first approach to this problem was developed among others by Andjel 
and Vares [2], who proved the hydrodynamic limit (1) for the asymmetric 
exclusion process when the initial datum is a single step with density A to 
the left and p to the right. This is usually referred to as the Riemann prob- 
lem for (1), in which case the solution has a simple explicit form when G 
is strictly convex or concave. Their approach combined this simplifying fea- 
ture and comparison with equilibrium systems, which are explicitly known. 
Comparison is allowed by the monotonicity (also called attractiveness) of 
the system. The method of Andjel and Vares [2] also applies to other attrac- 
tive systems with explicit invariant measures, such as the zero-range process; 
however in this case, one must add the assumption of flux concavity or con- 
vexity, which is not systematically true, nor easy to verify. Note that the 
fc-step exclusion process introduced by Guiol [12] is an attractive process 
with product invariant measures but simple nonconcave and nonconvex flux 
function. 

A breakthrough was achieved by Rezakhanlou [20] , who proved the hydro- 
dynamic limit for the asymmetric exclusion process in any dimension, with 
an arbitrary initial datum in (1). His approach is more general and based 
on entropy inequalities. It also applies, without any convexity assumptions, 
to other attractive processes with product invariant measures. 

In [5], we initiated a "resurrection" of the approach of Andjel and Vares 
[2]. Our goal was to show that, with a little more work, the same result 
could be reached in dimension one with the same approach as in [20], but in 
a constructive way. Our argument was based on two ideas. First, we relaxed 
the convexity assumption of Andjel and Vares [2] by introducing a variational 
formula for the Riemann solution. Next, we provided a general argument, for 
finite-range attractive processes, showing that the hydrodynamic limit for 
Riemann initial condition implied hydrodynamics for general initial datum. 
This argument is inspired by the Glimm's scheme in the theory of hyperbolic 
conservation laws (see, e.g., Chapter 5 of [25]), a procedure to reconstruct 
general entropy solutions from Riemann entropy solutions. 

In the present paper we refine our approach to obtain new results. Namely 
we establish the hydrodynamic limit for quite general finite-range attractive 
processes whose invariant measures are not product, nor even explicit at 
all, and not much is known about their properties: for instance the support 
of the densities of extremal elements is unknown. This extends two former 
results in this direction: the first one by Seppalainen [23] for the totally 
asymmetric, nearest-neighbor i'f -exclusion process, based on the author's 
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variational coupling method; the second one by Rezakhanlou [21] for the par- 
tially asymmetric nearest-neighbor iiT-exclusion process, using an abstract 
characterization of Hamilton- J acobi semigroups, and the ergodic theorem 
of Ekhaus and Gray [9] for nearest-neighbor attractive systems. 

In a forthcoming paper, we will demonstrate how a further refinement 
of our method leads to a strong law of large numbers (i.e., almost sure 
hydrodynamics) for attractive systems on Z. 

2. Notation and results. Throughout this paper N= {1,2,...} will de- 
note the set of natural numbers, Z + = {0,1,2,...} the set of nonnegative 
integers, and ffi + * = M + — {0} the set of positive real numbers. We consider 
particle systems on Z with at most K particles per site, K £ N. Thus the 
state space of the process is X = {0, 1, . . . , K} z , which we endow with the 
product topology. A function defined on X is called local if it depends on 
the variable rj £ X only through (r](x),x € A) for some finite subset A of Z d . 
If 7] is an X- valued random variable and v a probability measure on X, we 
write rj ~ v to specify that rj has distribution v. The notation v(f), where 
/ is a real- valued function and v a probability measure on X, will be an 
alternative for fxfdv. We write v n =? v to denote weak convergence of a 
sequence (i / ni n £ N) of probability measures on X to some probability mea- 
sure v on X, that is, f n (/) — > v{f) as n — > oo for every continuous function 
/onX. 

2.1. The model. The dynamics consists of particles' jumps, according to 
the Markov generator 

(2) Lf( v )= J2 p(y-v)Kv(x)My))lf(ri x ' y )-f(ri)] 

for a local function /, where rj x ' y denotes the new state after a particle has 
jumped from x to y [i.e., r] x ' y (x) = r/(x) — 1, r] x ' y (y) = rj{y) + 1, r] x ' y (z) = 
r](z) otherwise], p is the particles' jump kernel, that is, J2 Z £zP( z ) = 1j an d 
b : Z + x Z + — > M + is the jump rate. We assume that p and b satisfy: 

(Al) The semigroup of Z generated by the support of p is Z itself (stan- 
dard irreducibility); 

(A2) p is finite range, that is, there exists M > such that p(x) = for 
all \x\ > M; 

(A3) b(0,-) =0, b(-,K) = (no more than K particles per site), and 
b(i, j) > for < i < K and < j < K (nondegenerate jump rates); 
(A4) b is nondecreasing (nonincreasing) in its first (second) argument. 

Remark. In view of assumption (A4), the third condition in (A3) can 
be equivalently replaced by the simpler condition b(l,K — 1) > 0. 
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If b satisfies additional algebraic relations (see [8]), which include, for 
example, the asymmetric exclusion process, the system has an explicit family 
of invariant product measures indexed by the mean density of particles, and 
is the so-called misanthrope 7 s process. Here we do not assume these relations, 
and the system in general has no explicit invariant measures. This is even 
the case for the simplest possible b, that is, b(n,m) = l{ n >o, m <A'}> f° r which 
the system is called K-exclusion process — see [23] — (except for symmetric 
p, in which case there are still product invariant measures, see [13]; but then 
the relevant time scale is diffusive and not Euler, and Theorem 2.2 below 
yields a trivial nonevolving hydrodynamic limit). 

The special form 

(3) p(y-x)b(r](x),r)(y)) 

in (2), in which a jump kernel p(-) and a configuration-dependent part &(•, •) 
are decoupled, does not play any role in our paper, but we retained it by 
analogy with [8] , where it was needed to exhibit product invariant measures 
for a certain class of b's. Here we might as well consider more general jump 
rates of the form 

(4) b(y-x,r){x),r)(y)) 

for which the set of assumptions (Al)-(A4) would have to be replaced by: 

(Al') The semigroup of Z generated by the set {z£Z: inf n> o, m <_ft r b(z, n, 
m) > 0} is Z itself. 

(A2') Finite range assumption: there exists M > such that b(z, •, •) = 
for all \z\ > M. 

(A3') &(•, 0, •) = 0, &(•, -,K) = (no more than K particles per site). 
(A4') b is nondecreasing (nonincreasing) in its second (third) argument. 

Note that, in this more general context, assumption (Al') above replaces 
both assumption (Al) and the second part of assumption (A3). Under these 
assumptions, the coupling arguments of Cocozza-Thivent [8] needed in Sec- 
tion 3 below carry over. There is another reason why (3) is not so relevant 
here: for the misanthrope's process of Cocozza-Thivent [8], the decoupling 
in (3) is reflected in the hydrodynamic limit in the form of a decoupling in 
the macroscopic flux, that is 

(5) G(u)=jH(u), 

where 7 '■= J2 Z z p( z ) * s the mean drift depending only on p(-), and H(-) 
depends only on &(•). This is in fact due to existence of product invariant 
measures for all density values, see Remark 1. As a result, p(-) and &(•,•) 
have separate interpretations in the hydrodynamic limit. In our case, without 
product invariant measures, there is no a priori reason for a decoupling like 
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(5) in the macroscopic flux, hence no particular sense in the microscopic 
decoupling (3). 

Let the coordinatewise partial order on X be defined by n < £ if and 
only if rj(x) < £(x) for every x € Z. This induces a partial stochastic order 
for probability measures jjL\ and jj,2 on X; namely, we write [i\ < /i2 if the 
following equivalent conditions hold (see, e.g., [19]): 

(i) For every nondecreasing nonnegative function / on X, fj.\(f) < ^(f)- 

(ii) There exists a coupling measure ju on the product space X = X x X 
with marginals fj,\ and fj,2, such that ji{(r], £) : rj < £} = 1. 

It is shown in [8] that assumption (A4) implies the existence of a cou- 
pled Markov generator L (called the basic coupling) on X such that: (a) L 
generates Feller processes (?ft,£t)t>o whose components are Feller processes 
generated by L, and (b) for these coupled processes, r/o < £o &-s. implies 
Vt < £t a - s - f° r every t > 0. Properties (a) and (b) above imply monotonicity 
of the semigroup <S(-), that is, 

(6) ii<v => Vt>0,fj.S(t)<i>S(t). 

Either properties (a) and (b), or (6), are usually called attractiveness. These 
are still true for the more general model (4) under assumption (A4'). 

2.2. Scalar conservation laws and entropy solutions. We recall the defi- 
nition of entropy solutions to scalar conservation laws, which will appear as 
hydrodynamic limits of the above models. For more details, we refer to the 
textbooks of Godlewski and Raviart [11], Serre [25] or Bressan [6]. 

Let G : [0, K] — > M be a Lipschitz-continuous function, called the flux. It is 
a.e. differentiable, and its derivative G' is an (essentially) uniformly bounded 
function. We consider the scalar conservation law 

(7) d t u + d x [G(u)}=0, 

where u = u(x,t) is some [0, K] -valued density field defined on R x R + . 
Equation (7) has no strong solutions in general: even starting from a smooth 
Cauchy datum u(-,0) = uq, discontinuities (called shocks in this context) 
appear in finite time. Therefore it is necessary to consider weak solutions, 
but then uniqueness is lost for the Cauchy problem. To recover uniqueness, 
we need to define entropy solutions. 

Let (p : [0, K] — > R be a convex function. In the context of hyperbolic sys- 
tems, such a function is called an entropy. We define the associated entropy 
flux tp on [0,K] as 



i/j(u) := / 4>'{v)G'{v)dv, 
J o 
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((f>, ip) is called an entropy-flux pair. A Borel function u : R x R + * — ► [0, K\ is 
called an entropy solution to (7) if and only if it is entropy-dissipative, that 
is, 

(8) d t <p{u) + d x i/>(u) < 

in the sense of distributions on R x R + * for any entropy- flux pair (<p, ip). 
Note that, by taking cp(u) = ±u and hence if)(u) = ±G(u), we see that an 
entropy solution is indeed a weak solution to (7). This definition can be 
motivated by the following points: (i) when G and </> are continuously dif- 
ferentiable, (7) implies equality in strong sense in (8) (this follows from the 
chain rule for differentiation); (ii) this no longer holds in general if u is only 
a weak solution to (7); (iii) the inequality (8) can be seen as a macroscopic 
version of the second law of thermodynamics that selects physically rele- 
vant solutions. Indeed, one should think of the concave function h = — </> 
as a thermodynamic entropy, and spatial integration of (8) shows that the 
total thermodynamic entropy may not decrease during the evolution (this 
is rigorously true for periodic boundary conditions, in which case the total 
entropy is well defined). 

Kruzkov proved the following fundamental existence (Theorem 2 of [15]) 
and uniqueness (Theorem 5 of [15]) result: 

Theorem 2.1. Let uq:M — > [0,K] be a Borel measurable initial datum. 
Then there exists a unique (up to a Lebesgue-null subset o/Ix R + *j entropy 
solution u to (7) subject to the initial condition 

(9) lim u{-,t) =u o (0 in L} oc (R). 

This solution [has a representative in its L°°(R x R + *) equivalence class 
that] is continuous as a mapping 1 1— > u(-,t) from R + * to Lj oc (R). 

We recall here that a sequence (u n , n € N) of Borel measurable functions 
on R is said to converge to u in L 1 1 oc (R) if and only if 




for every bounded interval JcK. 

Remark. Kruzkov's theorems are stated for a continuously differen- 
tiable G. However the proof of the uniqueness result (Theorem 2 of [15]) 
uses only Lipschitz continuity. In the Lipschitz-continuous case, existence 
could be derived from Kruzkov's result by a flux approximation argument. 
However a different, self-contained (and constructive) proof of existence in 
this case can be found in Chapter 6 of [6]. 
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The uniqueness statement of Theorem 2.1 is in fact a consequence of the 
following L 1 -contraction and finite propagation properties (see also prop- 
erty (2) in Theorem 3.1 of [5], or [11] or [25]): 

Proposition 2.1. Let V := ||G'||oo denote the Lipschitz constant of the 
macroscopic flux G(-) in (7). Then, for any x < y and <t < (y — x)/{2V), 

ry—Vt ry 

(10) / \u(z,t) — v (z,t)\ dz < / \uq{z) — vq{z)\ dz, 

Jx+Vt Jx 

where u(-, •) and v(-,-), respectively, denote the entropy solutions to (7) with 
initial data uq(-) andvo(-). In particular, if Uq(-) and vq(-) coincide a.e. on 
[x,y], thenu(-,t) andv(-,t) coincide a.e. on [x + Vt, y — Vt] . 

We next recall a possibly more familiar definition of entropy solutions 
based on shock admissibility conditions, but valid only for solutions with 
bounded variation. This point of view selects the relevant weak solutions by 
specifying what kind of discontinuities are permitted. The following geomet- 
ric condition is known as Olemik's entropy condition (see, e.g., [11] or [25]). 
A discontinuity (u~ ,u + ), with := u(x ±0,t), is called an entropy shock, 
if and only if: 

The chord of the graph of G between u~ and u + lies: 

(11) below the graph if u~ < n + , above the graph if u + < 
u~~ . 

In the above condition, "below" or "above" is meant in wide sense, that is 
does not exclude that the graph and chord coincide at some points between 
u~ and u + . In particular, when G is strictly convex (resp. concave), one 
recovers the fact that only (and all) decreasing (resp. increasing) jumps are 
admitted. Note that, if the graph of G is linear on some nontrivial interval, 
condition (11) implies that any increasing or decreasing jump within this 
interval is an entropy shock. 

Condition (11) can be used to select entropy solutions among weak so- 
lutions. Let TV j denote the variation of a function defined on some closed 
interval / = [a, b] C K, that is, 

n-l 

TV/[n(-)]= sup ^2\u(x i+1 ) - u(xi)\. 

x =a<x 1 <---<x n =b j_Q 

Let us say that u = u(-,-) defined on R x M + * has locally bounded space 
variation if 

(12) supTVj[u(-,t)] < +oo 

teJ 

for every bounded closed space interval I C M and bounded time interval 
J C M + *. Then the following result, which will be used in Section 4.1, is a 
consequence of Vol'pert [26]. 
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Proposition 2.2. Let u be a weak solution to (7) with locally bounded 
space variation. Then u is an entropy solution to (7) if and only if, for a.e. 
t>0, all discontinuities of u(-,t) are entropy shocks. 

For completeness a proof of this statement is given in Appendix A. 

Remark. One can show that, if the Cauchy datum uq has locally bounded 
variation, the unique entropy solution given by Theorem 2.1 has locally 
bounded space variation. Hence Proposition 2.2 extends into an existence 
and uniqueness theorem within functions of locally bounded space variation, 
where entropy solutions may be defined as weak solutions satisfying (11), 
without reference to (8). 

2.3. The hydrodynamic limit. Before stating our main result, we recall 
some standard definitions in hydrodynamic limits. The integer N € N is the 
scaling parameter for the hydrodynamic limit, that is, the inverse of the 
macroscopic distance between two consecutive sites. The empirical measure 
of a configuration rj viewed on scale N is given by 

(13) a N ( V ,dx) = N" 1 £ V(y)$y/N(dx) G M, 

where M denotes the set of positive, locally finite measures on R. The set A4 
is equipped with the topology of vague convergence, defined by convergence 
for continuous test functions with compact support. Let u(-) be a determin- 
istic bounded Borel function on R. A sequence of random configurations r) N 
is said to have density profile u(-) if a N (r] N ,dx) converges in probability to 
u(-) dx as N — > oo. 

Let u(v) be a deterministic bounded Borel function on R x R + . A se- 
quence of processes r/^ = (77^ , t > 0) generated by L, with random initial 
configurations t]q , is said to have hydrodynamic limit u(-,-) (under Euler 
time scaling), if for every t > 0, the sequence has density profile u(-,t). 

We now state our main result. 

Theorem 2.2. Assume the sequence t]q has density profile uq(-), where 
uo(-) is a measurable [0, K] -valued profile on R. Then the sequence of pro- 
cesses rj N has a hydrodynamic limit u(-,-) given by the unique entropy so- 
lution to the scalar conservation law (7) with initial condition uq, where 
G : [0, K] — > R + is a Lipschitz- continuous flux function defined below. 

The flux function G in Theorem 2.2 is obtained as follows. Let J and 
iS denote respectively the set of invariant measures for L, and the set of 
shift-invariant measures on X. Since X is compact, and the process has the 
Feller property, J n S is a nonempty (the empty and full configurations are 
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unchanged by the dynamics) convex compact subset of the set of probability 
measures on X. Thus any element v € J n S is a mixture of its extremal 
elements. In the usual setting of Andjel and Vares [2], Rezakhanlou [20], 
Bahadoran et al. [5], it is known that 

(J n«S) e = { W £[(),#]}, 

where v p are product measures and the index e denotes the extremal ele- 
ments. In our setting, there are no explicit invariant measures. However, as 
we will explain in Section 3, the same coupling arguments as used by Liggett 
[18], relying only on attractiveness and irreducibility, are still sufficient to 
establish 

(14) (3ns) e = {u p , P en}, 

where TZ is a closed subset of [0,K] containing and K, and v p is a shift- 
invariant probability measure on X with fp[f?(0)] = p. Whether TZ = [0,K] 
is an open problem. 

We define the microscopic flux (here across site 0) by 

(15) j{rj)= p(y- x ) b (v( x ),v(y)) - p(v - x ) b (v( x ),v(y))- 

x<0,y>0 x>0,y<0 

The term "microscopic flux" is justified by the microscopic conservation 
equation [clear from (2)] 



(16) L 



r x j(rj) -T y j(ri) 



for x < y, where r is the space shift, that is, T x r\(z) = r\(x + z) and T x j(rj) = 
j(r x r)). Equation (16) says that the instantaneous variation of the number 
of particles in interval (x,y] n Z is due to the fluxes (translates of j) at the 
boundaries. The macroscopic flux is now defined by: for p € TZ, 

(17) G(p) = u p [j(v)], 

then interpolate G linearly on the complement of TZ, which is an at most 
countable union of disjoint open intervals. Note that, for p 6 TZ, by shift- 
invariance of v p , we also have G(p) = v P \j'{rj)\ with j'(rj) = J2 Z z P{ z )^{ r l{^) ^V^)) ■ 



Remark 1. For the misanthrope's process, the decoupled form (5) fol- 
lows from (15), (17), and the fact that v p exists and is a product measure 
for all p. Indeed, we can then define H(p) := u p [b(r](x),r](y))] for any x ^ y, 
because this quantity is independent of x and y. Hence (5) holds with 

7 : = p{y-x)- Y p(y- x )=Y z P( z )- 

x<0,y>0 x>0,y<0 z 
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Remark 2. Unbounded systems. The arguments and results of Section 3 
do not use the assumption that the number of particles per site is bounded by 
K < +oo. Thus they are still true for systems with state space N z satisfying 
assumptions (Al), (A2), (A4) and 

(A3") 6(0, •) = and b(i,j) > for % > 0. 
In this case (14) still holds with 1Z a closed subset of M + . However, in order 
to define a macroscopic flux by interpolation, we would need to know that 7Z 
is unbounded, which we cannot establish. Then the arguments and results 
of Section 4 would also extend, thus establishing Riemann hydrodynamics. 
On the other hand, we do not currently know how to carry out the passage 
from Riemann to general hydrodynamics (Section 5) for unbounded systems 
without product invariant measures, because we cannot prove Lemma 5.1 
for such systems. 

Remark 3. Nothing can be said in general about the set 1Z, except 
that it is closed and contains and K. However for the totally asymmetric 
nearest-neighbor K-exclusion process, one can benefit from some properties 
of G derived in [23], to get slightly more precise information about 1Z. 

Corollary 2.1. For the totally asymmetric K -exclusion process, and 
K are limit points oflZ, and 1Z contains at least one point in [1/3, K — 1/3]. 

Proof. The following properties were established in [23] for G: 
(a) G is symmetric around u = K/2: G(K — u) = G{u) for every u 6 



(b) G is concave. 

(c) G has the following bounds: F(u) < G(u) < H (u) for every u £ [0, K], 



If were no limit point for 1Z, this would imply by construction of G that 
there exists e > such that G is linear on [0,e]. But there can be no linear 
function lying between the lower bound F and upper bound H near the 
origin, because F and H have the same derivative at u = and are strictly 
concave near the origin. The same argument works around K by symmetry 
of F and H. 

Let us now assume 1Z n [1/3, K — 1/3] = 0. Since 7Z is closed, we must 
have K n [1/3 - e,K - 1/3 + e] = for some < e < 1/3. Then, by con- 
struction of G, the graph of G between u = 1/3 — e and u = K — 1/3 + e 



[0,K]. 



where 




if 0<u<l/2, 

if l/2<u<K -1/2, 

if K- 1/2 <u<K, 




u) 



if < u < K/2, 
if K/2 <u<K. 
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must be a straight line. By the symmetry property (a) of G, this line must 
be horizontal, that is, G has constant value over [1/3 — s, K — 1/3 + s]. Be- 
cause of the lower bound in (c), this constant value must be at least 1/4. 
Hence G(l/3 — e) > 1/4, which contradicts the upper bound in (c), since 
#(1/3 -e) < #(1/3) = 1/4. □ 

2.4. Riemann solutions and Glimm's scheme. Of special importance 
among entropy solutions are the solutions of the Riemann problem, that 
is, the Cauchy problem for particular initial data of the form 

(18) uq(x) = Xl{ x<0 } + pl{x>o}- 

Indeed: (i) these solutions can be computed explicitly and have a variational 
representation: see Section 4.1; (ii) one can construct approximations to the 
solution of the general Cauchy problem by using only Riemann solutions. 
This has inspired our belief that one could derive general hydrodynamics 
from Riemann hydrodynamics. 

We will briefly explain here the principle of approximation schemes based 
on (ii), the most important of which is probably Glimm's scheme, introduced 
in [10]. Consider as initial datum a piecewise constant profile with finitely 
many jumps. The key observation is that, for small enough times, this can 
be viewed as a succession of noninteracting Riemann problems. To formalize 
this, we recall part of Lemma 3.4 of [5], whose proof was a direct application 
of Proposition 2.1 above. We denote by R\ tP (x,t) the entropy solution to 
the Riemann problem with initial datum (18). 

Lemma 2.1. Let xq = —oo < x\ < • ■ • < x n < x n+ \ = +oo, and e := 
minfe(xfc + i — Xk). Consider the Cauchy datum 

n 

u :=$> fe l (a;fcja;fc+l) , 

k=0 

where r^ € [0, #]. Then for t < e/(2V), with V given in Proposition 2.1, 
the entropy solution u(-,t) at time t coincides with R rk _ lt r k (' ~ x k,t) on 
(xk—i + Vt, Xk+i — Vt). In particular, u(-,t) has constant value r& on (x& + 
Vt,x k+1 -Vt). 

Given some Cauchy datum uq, we construct an approximate solution 
tt(-,-) for the corresponding entropy solution u(-,-). To this end we define 
an approximation scheme based on a time discretization step At > and a 
space discretization step Ax > 0. In the limit we let Ax — > with the ratio 
R := At/Ax kept constant under the condition 



(19) 



R<1/(2V) 
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known as the Courant-Friedrichs-Lewy (CFL) condition. Let tk '■= kAt de- 
note discretization times. We start with k = 0, setting Uq :=Uq. 

Step 1 (approximation step): Approximate with a piecewise constant 
profile ut whose step lengths are bounded below by Ax. 

Step 2 (evolution step): For t € [£fc,tfc+i), denote by %(■,£) the entropy 
solution at time t with initial datum at time tk- By (19) and Lemma 2.1, 
Uk(-,t) can be computed solving onlyRiemann problems. Set u^ +1 = ife+i). 

5"tep 3 (iteration): Increment fc and go back to step 1. 

The approximate entropy solution is then defined by 

(20) u(-,t) :=J2u k (;t)lit k ,t k+1 ){t). 

The efficiency of the scheme depends on how the approximation step is 
performed. In Glimm's scheme, the approximation ut is defined as 

(21) u+ := %((j + a fc/2)Ax)l ((j _ 1/2 )A Xi ( j+ i/2)A2), 

jek/2+z 

where S (—1,1). Then we have the following convergence result. 

Theorem 2.3. Let uq be a given measurable initial datum. Then ev- 
ery sequence e n [ as n — > oo has a subsequence S n | such that, for a.e. 
sequence (a^) w.r.t. product uniform measure on (— 1, l) z , the Glimm ap- 
proximation defined by (2Q) and (21) converges to u in Lj oc (R x M + *) as 
Ax = 5 n i 0. 

When no has locally bounded variation, the above result is a special- 
ization to scalar conservation laws of a more general result for systems of 
conservation laws: see Theorems 5.2.1, 5.2.2, 5.4.1 and comments following 
Theorem 5.2.2 in [25]. In Appendix B, we prove that it is enough to assume 
uq, measurable. 

Due to the nature of the approximation step (21), the proof of Theorem 2.3 
does not proceed by direct estimation of the error between v/^ and u(-,tk), 
but indirectly, by showing that limits of the scheme satisfy (8). Our approach 
is based on a different approximation, introduced first in Lemma 3.6 of [5], 
and refined in Lemma 5.5 of the present paper. This approximation allows 
direct control of the error by using a suitable distance, defined and denoted 
by A(u,v) in Section 5. Intuitively errors accumulate during approximation 
steps, but might be amplified by the evolution steps. The key properties 
of our approximation are that (i) the total error accumulated during the 
approximation step is negligible as e — > 0; (ii) the error is not amplified by 
the evolution step, because the resolution semigroup of (7) is A-contractant; 
(iii) a similar property holds at particle level, that is, the particle system 
is contractant for a microscopic version of A. This allows us to mimic the 
scheme at particle level. 
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3. Properties of invariant measures. The main result of this section is 

Proposition 3.1. (i) (14) holds with 1Z a closed subset of [0,K] con- 
taining and K , and v p a stationary shift-invariant measure such that 

(ii) The measures v p are stochastically ordered: v p < v p < if p< p' . 
(hi) v p has a.s. density p, that is, (21 + 1) _1 Ylx=-i v( x ) converges v p -a.s. 
to p as I — > oo. 

Remark. A proof of (i) and (ii) can be found in [9] and [21] for nearest- 
neighbor processes. It is extended [as well as (hi)] to finite-range if -exclusion 
processes in [24]. In these papers it appears as a byproduct of a large time 
convergence result, which also implies an ergodic theorem. However if one is 
only interested in (i) and (ii), and not in the ergodic theorem, a much easier 
independent proof is possible, along the lines of [18] or [19], without using 
any convergence result. This is what we explain below. 

If the elements of (In5) e were product measures, it would immediately 
follow that they are ergodic with respect to translations, implying (hi). Since 
this is not the case here, for a proof of (hi) we refer the reader to [21]; 
though in that paper the proof is given in the context of nearest-neighbor 
if -exclusion process, the same argument shows that, for any conservative 
finite-range system, every element of (Ifl5) e has an a.s. density. 

We divide the proof of (i)-(ii) into two lemmas. 

Lemma 3.1. Any two elements p>\ and p2 of (IC\S) e are ordered stochas- 
tically as their mean densities, that is: p\ < pi if p\\r\(fS)\ < p2[v(®)]- ^ n par- 
ticular, for every density p G [0,if], there is at most one v G (TnS) e with 
v[ri(0)] = p; when it exists, we denote it by v p . 

Proof. The proof uses standard coupling arguments, which we briefly 
recall for self-containedness. We use couplings on the space X and denote 
with a "tilde" quantities related to the coupled process. For any p\ and 
P2, elements of (In5) e , there exists a coupling measure p with marginals 
pi and p2 which belongs to the set (Z n S) e of invariant and translation 
invariant extremal measures for the coupled process (see [18] and [19]). By 
Theorem 1.9 in [8], p is supported on ordered configurations, that is, (r/,£) 
such that r] < £ or £ < rj. 

Remark. To be precise we recall (as mentioned above in Section 2.1) 
that the difference between our model and the misanthrope's process studied 
in [8] is that for the latter additional assumptions are imposed on b to ensure 



14 



BAHADORAN, GUIOL, RAVISHANKAR AND SAADA 



existence of product invariant measures. However these assumptions are not 
required for the proof of Theorem 1.9 in [8], which only uses assumptions 
(Al), (A3) and (A4), and extends in a natural way the coupling arguments 
introduced by Liggett [18] for the simple exclusion process. For details of 
the extension to misanthrope's process see Section 5 of [1] and [7]. 

It follows from extremality of p that the event {?/ < £} has ^-probability 
or 1 (cf. Chapter 8, Proposition 2.13 of [19]). This and pi[r)(0)] < /i 2 [f (0)] 
imply #({77(0) < £(0)}) = 1. Indeed, £({77(0) > £(0)}) > would imply £({77 < 
£}) = 0, hence //({£ < 77}) = 1; but the latter combined with //({77(G) > 
£(0)}) > in turn implies the contradiction px[r)(0)] = /x[t?(0)] > /2[£(0)] = 
A*2[£(0)]. Now, since p, is translation invariant, p({r)(Q) < £(0)}) = 1 implies 
Jl({i](x) < £(cc)}) = 1 for all x G Z. Hence p({rj < £}) = 1, and thus p\ < pi- 
□ 



Lemma 3.2. TZ is a closed subset of [0, K]. 



Proof. Since the empty and full configurations are unchanged by the 
dynamics, TZ contains at least and K. Let (p n )n>o be an increasing se- 
quence in TZ, converging to p (with obvious changes, what follows can be 
adapted to a decreasing subsequence). By Lemma 3.1, (f Pn ) n >o is an in- 
creasing sequence in (IC\S) e , therefore converges (weakly) to some u; it is 
easy to see that v belongs to X n S. Suppose p ^ TZ, so v is not extremal. 
Then v = Jj K , v a r y(da) for some measure 7 on [0,K] such that 7({p}) = 0, 
and 



v[ri(0)]= lim v [rj(0)]=p= / aj(da). 

rwoo J[0,K] 

Then, there exists some p' < p such that "f([0,p']) > 0, and by (iii) of Propo- 
sition 3.1, 



(22) vlr, GX: lim (2Z + l)" 1 V rj{ 

\ x=—l 



a?) <p'j >0. 

Let n be such that p n > p' , then because u Pn < u, there exists a coupling 
measure z7, with marginals v and u„ n , such that ?({(t?,£) G X:£ < 77}) = 1. 
This implies [using (iii) of Proposition 3.1] 

1 1 
lim (2/ + l)" 1 V n(x) > lim (21 + l)" 1 V £(x) = p n , x7-a.s. 

i— >oo ' — ', !— >oo z — ' 

x=—l x=—l 

Therefore v({rj G X: lirm_ >00 (2i + l)" 1 Z)L=-j ^O^) — Pn}) = 1, which contra- 
dicts (22). Thus v is extremal with density p. So p belongs to TZ, proving 
that TZ is closed. □ 
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4. The 7?--valued Riemann problem. In this section we derive local equi- 
librium for the so-called Riemann problem. This is a first step to obtain 
Theorem 2.2, for the particular initial condition (18) with A,p € 7£. In other 
words we generalize [5], Section 2. We first define a corresponding initial 
measure v\ %p for the particle system which, as opposed to Andjel and Vares 
[2] and Bahadoran et al. [5], is generally not product. The measure u\ iP is 
defined in such a way as to enjoy the following properties: 

(i) Negative (nonnegative) sites are distributed as under v\ {v p )\ 

(ii) Tiv x , P > v\, P (tiV\, p < v\ p ) if A < p (A > p); 

(iii) v\ )P is stochastically increasing with respect to A and p. 

In [5], v\ sP was product, and given explicitly by its marginals site per site. 
In the present setting we construct a (nonexplicit) v\ p using the fact that 
if A < p, then v\<v p . Let Jibe a coupling measure on X for v\ and v p such 
that /x({(f7,£) : f] < £}) = 1. We define the measure u\ sP as the distribution of 
the X-valued random variable 7 = 7(77, £) defined on the coupling probability 
space (X, p) by j(x) = rj(x) if x < and j(x) = £(x) if x > 0. 

We shall prove in Section 4.2 below that the system starting from v\„ 
is in local equilibrium with density profile given by the entropy solution to 
the Riemann problem with initial datum uo(-). As a consequence we will 
obtain, in Corollary 4.2, the hydrodynamic limit for the 7£-valued Riemann 
problem. As a first step, we prove in Section 4.1 that the 1Z- valued Riemann 
problem still admits the same variational characterization as in [5]. 

4.1. A variational principle for Lipschitz-continuous flux. The minimum 
regularity we must ensure for G in order to extend the variational charac- 
terization of Bahadoran et al. [5] is Lipschitz continuity. We therefore prove 
the following: 

Lemma 4.1. The flux G defined in (17) is Lipschitz-continuous. 

Proof. It is easy to see that there is a constant C > and a finite set 
5cZ such that 

m-m\<cj2\v(x)-c(x)\. 

For p, r € 1Z, with p < r, we have v p <v T , so there is a coupling measure v p ^ r 
for v p and u r under which 77 < £ a.s.; thus 

\G{p) - G(r)\ < v p , r (\j(v) - J(C)I) < C\S\(r - p). 

Hence G is uniformly Lipschitz inside 7Z, thus everywhere, since it is com- 
pleted by linear interpolation and 1Z is closed. □ 
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From now on we assume A < p; adapting to the case A > p is straightfor- 
ward, by replacing in the sequel lower with upper convex hulls, and min- 
ima/minimizers with maxima/maximizers. Consider G c , the lower convex 
envelope of G on [X,p]. There exists a nondecreasing function H c (hence 
with left /right limits) such that G c has left /right hand derivative H c (a±0) 
at every a; the notation H c (a ± 0) means left /right hand limit of H c at a. 
The function H c is defined uniquely outside the at most countable set of 
nondifferentiability points of G c 

(23) G = {a G [A, p] : H c {a - 0) < H c (a + 0)}. 

As will appear below, the particular choice of H c on does not matter. 
Let = t>*(A,p) = H c {\ + 0) and v* = v*(X,p) = H c (p — 0). Since H c is 
nondecreasing, there is a nondecreasing function h c on [i>*,w*] such that, for 
every v G [v*,v*], 

a < h c (v) H c (a)<v, 

(24) 

a>h c (v) H c (a)>v. 

Any such h c satisfies 

h c (v - 0) = inf {a € R : H c (a) >v} = sup{a G R : H c {a) < v}, 

(25) 

h c {y + 0) = inf {a G M : H c (a) > v} = sup{a G R : H c (a) < v}. 

It is not difficult to see that, anywhere in (25), H c {a) may be replaced with 
H c (a ± 0). The following properties can be derived easily from (24) and (25): 

1. Given G, h c is defined uniquely, and is continuous, outside the at most 
countable set 

Siow(G) = {v G [t>*,v*] :G C is differentiable with derivative v 

(26) 

in a nonempty open subinterval of [A,p]}. 

By "defined uniquely" we mean that for such v's, there is a unique h c (v) 
satisfying (24), which does not depend on the choice of H c on 0. 

2. Given G, h c (v ± 0) is uniquely defined, that is, independent of the choice 
of H c on 0, for any v G For v G Si ow (G), (h c (v — 0), h c (v + 0)) is 
the maximal open interval over which H c has constant value v. 

3. For every a G and v G (H c (a — 0),H c (a + 0)), h c (v) is uniquely defined 
and equal to a. 

In the sequel we extend h c outside [t>*,i>*] in a natural way by setting 
(27) h c (v ) = A fovv<v*, h c (v) = p (or v > v* . 

Proposition 4.1. Let A, p G TZ, A < p. 
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(i) For every v G M\Si ow (G), G(-) — v- achieves its minimum over [A, p] 
at the unique point h c (v), and u(x,t) := h c (x/t) is the weak entropy solution 
to (7) with Riemann initial condition (18). 

(ii) The previous minimum is unchanged if restricted to [A,p] CiTZ. As a 
result, the Riemann entropy solution is a.e. TZ-valued. 

This extends Proposition 2.1 of [5], where we assumed G € C 2 (M). In the 
present setting part of the proof must be modified, because it explicitly used 
the regularity assumption. However the following lemma (Lemma 2.2 of [5]) 
carries over with the same proof. 

Lemma 4.2. For w£l, G(-) — v ■ and G c {-) — v- have the same global 
minimum value on \\,p\; the set of global minimizers of G c {-) — v- on [X, p] 
is the interval [h c (v — 0), h c {v + 0)]; h c (v — 0) and h c (v + 0) are, respectively, 
the smallest and greatest global minimizer of G(-) — v-. 

As a consequence of Lemma 4.2 we obtain the following: 

Corollary 4.1. (i) If v e T>\ ow (G), the graph of G c between h c (v - 0) 
and h c (v + 0) is a chord of the graph ofG, with slope v, lying below the graph 
ofG. 

(ii) The graphs of G and G c coincide at points a = h c (v), where v £ 
R\Eiow(G)- 

(hi) Assume a€ [A,p] and G c (a) =G(a). Then the restriction to [A, a] 
(resp. [a,p]) of G c is the lower convex envelope of the restriction of G to 
[A, a] (resp. [a,p]). 

PROOF, (i) Let v e T,\ ow (G). By Lemma 4.2, the points (h c (v ± 0), 
G c (h c (v ± 0))) lie on the graphs of G and G c - By convexity, G c (-) — v- is 
constant between the two global minima h c (v ±0), so the graph of G c be- 
tween these points is linear with slope v. It is thus also a chord for G, and 
lies below G because G C <G. 

(ii) Follows immediately from continuity of h c at v and Lemma 4.2. 

(hi) Denote by G\ (resp. G^.) the lower convex hull of G on [A, a] (resp. 
[a,/?]). Then G c < G\ on [A, a], since on this interval G c is convex and 
G c < G. Similarly we have G c < G 2 C on [a,p]. We now prove the converse in- 
equalities. At endpoints the lower convex hull of a continuous function must 
coincide with the original graph (otherwise one could increase the convex hull 
in the neighborhood of an endpoint where it does not coincide, while keep- 
ing it convex and below G). Hence G\{a) = G 2 c {a) = G(a) = G c (a). Denote 
by H l c (a ± 0) the left /right hand derivative of G\ at a for z = 1,2. Because 
(as argued above) G\ > G c on [A, a], we must have (a — 0) < H c (a — 0). 
Similarly H 2 C (a + 0) > H c (a + 0) , and thus Hi (a - 0) < H 2 (a + 0) . It follows 
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that the function G cc obtained by concatenation of G\ and G\ is convex on 
[A, p]. Since G cc < G, we must have G cc < G c , which implies G\ < G c (resp. 
G 2 C < G c ) on [A, a] (resp. [a,p]). 

□ 



Proof of Proposition 4.1. (i) We set u{x,t) :=h c (x/t) for x/t ^ 
Si ow (G). Since h c is nondecreasing, it has bounded variation; it follows that 
u has locally bounded space variation in the sense of (12). To prove that u 
is the unique entropy solution given by Theorem 2.1 for Cauchy datum uo, 
we check points (a), (b), (c) below: indeed the first two of them, combined 
with Proposition 2.2, imply that u is an entropy solution: 

(a) u is a weak solution, 

(b) jumps satisfy Olemik's entropy condition (11), 

(c) u(-,t) has initial datum uq in the sense of (9). 



Proof of (a). Taking test functions, we see that u(x,t) = h c (x/t) is a 
weak solution of the hydrodynamic equation (7) if and only if h c (-) satisfies 

(28) ^-[G(h c (v))-vh c (v)] = -h c (v) 

av 

weakly with respect to v. From Lemma 4.2, 

G(h c (v)) - vh c (v) = inf [G(-) - v] = inf [G c (-) - v] = -G* c (v) 
[A,p] [A,jo] 

for a.e. v € M, where G* is the Legendre-Fenchel transform of G c (see, e.g., 
[22]), therefore also a convex function, and (G*)' = h c (v ) almost surely with 
respect to Lebesgue measure. Using the absolute continuity of G* (which 
follows from convexity) (28) is obtained through integration by parts. □ 



Proof of (b). This follows from (i) of Corollary 4.1. □ 



Proof of (c). By (27), u(-,t) coincides with u outside [u*t, v*t]. □ 



(ii) Assume that, for some v ^ T,\ ov/ (G), the unique minimum of G(-) — v- 
is achieved at a ^ 1Z. Then a € (A,/)) and there exists an open interval 
/ C (A,p) containing a such that G is linear with slope w ^ v on / (recall 
that G is linearly interpolated on 1Z C ). But then G(-) — v- cannot achieve a 
minimum inside /. □ 



4.2. Local equilibrium. We now state the main result of this section: 
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Proposition 4.2. For every v £l\ £i ow (G7), as t — > oo, 

where u(x,t) = u(x/t, 1) is the entropy solution to (7) constructed in Propo- 
sition 4.1 for the initial datum (18). 

A consequence is the following: 

Corollary 4.2. Assume t]q ~ v\ p for every N 6 N. Then the sequence 
of processes rj^ = (n^ , t>0) has hydrodynamic limit u(-, •) as above. 

Remark. This is not an application of Landim [16], since we are no 
longer dealing with product measures. 

Proof of Corollary 4.2. Pick a continuous test function / with 
compact support in R, and I £ N. Then 



f f(x)a N ( V , dx) = N- 1 f(x/N) V ( 



N- 1 Y,f(x/N)r ] l (x)+o N (l) 



f(x)rj l ([Nx])dx + o N (l) 



where rf(x) = (21 + 1) 1 J2\ y - X \<l v(y)i an d 0^(1) is a vanishing sequence 
depending on / but not on n. Hence, 



(29) 

< J f{x)v Kp T [Nx] S{Nt)(\rf{U)-u{x/t,l)\)dx + o N {l) 
By Proposition 4.2 and (hi) of Proposition 3.1, 



lim lim v x , p T [Nx] S(Nt)(\rj l (0) - u(x/t,l)\) 

(30) 



= lim v u(x/t ,i)(W(0) ~ u(x/t, 1)|) = 

for a.e. The result follows from (29), (30) and Lebesgue's theorem. 

□ 



The proof of Proposition 4.2 is based on a scheme introduced by Andjel 
and Vares [2] for one-dimensional attractive processes with strictly concave 
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(or convex) flux function G and for which (I D S) e consists of a contin- 
uous family of product measures indexed by their mean densities. In [5], 
we extended the proof up to nonconvex, nonconcave flux functions via a 
new variational formula. We generalize the latter proof, taking into account 
1Z. Therefore we sketch the arguments that do not need modifications, and 
detail the others. 

Proof of Proposition 4.2. We assume \,pe1Z and < A < p (the 
case A > p is similar). The first step consists in showing that the limit in 
Proposition 4.2 holds in Cesaro sense. As in [5], we can transpose Lemmas 
3.1 and 3.2 of [2] to our setting to get: 

Lemma 4.3. (i) Any sequence T n — > oo has a subsequence T Um for which 
there exists a dense set Del such that for any v G D 



where F(w) = f[G(a) — wa\ r ) w {da) for w G D. 

The proof of (i) uses attractiveness, (i) and (ii) of Proposition 3.1, and 
the fact that the measure v\ iP satisfies v\ < v\ jP < T\u\ tP < v p \ (ii) also uses 
the microscopic conservation equation (16). For these reasons the proofs of 
Andjel and Vares [2] carry over here without modification, as they made 
no use at that level of additional properties or assumptions (i.e., product 
invariant measures for all densities and strictly concave or convex flux G). 

Next we can show, as in [5]: 

Lemma 4.4. For all v £ T,\ ow (G), 



This is equation (9) from [5], except that here the minimum is restricted 
to 71. In fact we can exactly reproduce the proof given on pages 9-10 of 
[5]. Indeed, this proof uses only (a) attractiveness, (b) Lemma 4.3 above, 
(c) comparison of the process starting from u\^ p with the one starting from 
vq p for intermediate densities 9 G [A, p] , and (d) comparison of the process 




for a probability measure "f v on [X,p] C\1Z. 
(ii) Moreover for any u <v G D 

/ [«T» m ] 




(31) 
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starting from vq_ p with the one starting from uq. The only difference is that, 
in the present setting, we are constrained to choose 9 6 1Z; this is why the 
restriction to 7Z appears in (31). However from (ii) of Proposition 4.1, the 
infimum on the r.h.s. can be equivalently taken on the whole interval, and 
whenever it is uniquely achieved, this occurs in 1Z. Hence, 

lv = & u {v,i) VwGlR\S low (G) 

which concludes the first step. The second step is to replace the Cesaro limit 
by convergence in distribution. To this end, as in [5], we need Proposition 
3.5 from Andjel and Vares [2]. This proposition uses only attractiveness, so 
it can be repeated here. 

Lemma 4.5. Assume fx satisfies: 

(&) li < v\, (h) llt\ > li, (c) there exists vq finite so that 

1 f T 

hm - / LiT [vt] S{t)dt = u x 

for all v > vq . Then 

lim LiT\ vt ]S(t) = v\ for all v > vq. 

This immediately implies 

lim v\ p Tr vt ]S(t) = v\ for all v < v*, 

(32) 

^lim v\ :P T^ vt ]S(t) = v p for all v > v* . 

Now consider v £ [i>*,u*] \ Ei ow (G). Set a = h c (v). Let {v n ) n be a sequence 
such that v n < v, v n — > v and v n T,\ ow (G) [the latter point is possible since 
Siow(G) is countable]. Set a n = h c (v n ); then a n 6 TZ by (ii) of Proposi- 
tion 4.1. By monotonicity of h c , a n < a; by definition of h c , H c (a n — 0)<v n 
and thus H c (a n — 0) < v; and by continuity of h c at v, a n — > a. 

Let u n (x,t) =u n (x/t, 1) denote the entropy solution to (7) with initial 
datum Mq(x) = Al{ x<0 } + a n l{ 2 .> }, and G n (resp. G?™) the restriction of G 
(resp. G c ) to [A,a n ]. Since v n T,\ ow (G), by (ii), (hi) of Corollary 4.1, is 
the lower convex hull of G n . Hence, for every w > H c (a n — 0) = v*(X,a n ), 
we have, by (32), v\,a n T[ w t]S(t) — > v an weakly for every such w. In partic- 
ular we have v\. an T[ v t]S(t) — ► v an weakly. By attractiveness, v\ tP T\ vt ]S(t) > 
v\,a n T\vt]S(t). Therefore any subsequential weak limit ji v of i/\ tP T[ vt }S(t) sat- 
isfies jl v > v an ; since a n — ► a we get Ji v >v a - To get the reverse inequality 
consider a sequence v n > v, v n — > v with v n ^ £i ow (G7), set a n = h c (v n ) and 
consider □ 
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5. General hydrodynamics. In the previous section (Corollary 4.2) we 
established the hydrodynamic limit for certain initial measures associated 
with 1Z- valued Riemann initial profiles. In this section we prove that this 
implies Theorem 2.2, that is, the hydrodynamic limit for any initial sequence 
associated with any measurable initial density profile. To this end we adapt 
Section 3 of [5]. In that paper, we first proved Riemann hydrodynamics 
for any left- and right-hand densities. We then proceeded in the spirit of 
Glimm's scheme described in Section 2.4, by which one constructs general 
entropy solutions using only Riemann solutions. In our case we could only 
prove Riemann hydrodynamics for 1Z- valued left and right-hand densities, 
where 1Z is a closed subset of [0, K] outside which the macroscopic flux G in 
(7) is linear. Therefore we refine the procedure in such a way as to reconstruct 
a general entropy solution by using only 7£-valued Riemann solutions. 

To extend hydrodynamics from Riemann to general initial profiles, we 
pointed out in [5] that the main property required for the particle system 
was what we called macroscopic stability. To state this property (Lemma 5.1 
below) we need the following notation. Let u(-) and v(-) be two integrable 
[0, K]- valued density profiles on R, and 77, £ be two particle configurations 
in X with finitely many particles. Then we set 



A(«(-)>0) := sup 



(33) 



A N (r),£) = AT 1 sup 



A N (r],u(-)) := AT 1 sup 

x&L 



[u(y) - v(y)\ dy 

£fo(v)-«v)] 

y<x 



y<x 



x/N 



u{y) dy 



Lemma 5.1. Consider a sequence of coupled processes [t]^ ,C^) ( v ^ a the 
basic coupling) such that 



su P N- 1 ^2[ V ^(x) + ^(x)]<+o , 



a.s. 



Then, for every t > 0, 

(34) A N (r,% t ,C% t )<A N ( V »X?)+o N (l) 

where ojv(1) denotes a sequence of random variables converging to in prob- 
ability. 



A consequence of (34) is that the hydrodynamic limit depends only on 
the density profile at time 0, and not on the underlying microscopic struc- 
ture: that is, if a hydrodynamic limit holds for some initial sequence with 
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given density profile, it still holds for any initial sequence with the same 
density profile. For the system defined by (2), (34) follows from the argu- 
ments of Bramson and Mountford [4], Proposition 3.1 (the latter is written 
for the asymmetric exclusion process; however the only properties involved 
are attractiveness, irreducibility and limited number of particles per site, so 
the proof can be adapted to our model). In the nearest-neighbor case [i.e., 
when p(l) + p(— 1) = 1], (34) holds without even the on(1) and is only a 
consequence of attractiveness, that is, assumption (A4). It follows also from 
a slightly more general result, that is, Proposition 2.1 of [3]. 

Another property we shall need is the finite propagation property for the 
particle system (see also Lemma 3.1 of [5]), which follows from the finite- 
range assumption on p(-). This is a microscopic counterpart to the finite 
propagation property for entropy solutions (Proposition 2.1). 

Lemma 5.2. For any x,yEl*, any coupled process (r]t,Ct)t>o (via the 
basic coupling), and any < t < (y — x)/{2V'): if ijq and Co coincide on the 
site interval [x,y], then i]t and Ct coincide on the site interval [x + V't,y — 
V't] n 7L with probability at least 1 — e~ ct , where V and C are constants 
depending only on &(•,•) andp(-). 

The proof of Theorem 2.2 from Corollary 4.2 can be decomposed into two 
steps. The first step is to prove the result for 1Z- valued entropy solutions. 
This step follows essentially the proof of Theorem 3.2 from [5], with some 
refinements due to 1Z. The second step, specific to this paper, is to show 
that an arbitrary entropy solution can be approximated by 1Z- valued entropy 
solutions. In the following, we shall recall the main steps followed in [5], but 
mainly insist on the differences that arise in the present setting. 

5.1. Hydrodynamics for IZ-valued solutions. In Theorem 3.2 of [5], we de- 
rived general hydrodynamics from Riemann hydrodynamics. We can repeat 
the same arguments here, when restricted to 7£-valued entropy solutions. 
Before doing so, we need the following additional result: 

Lemma 5.3. Assume the initial datum uq{-) is a.e. IZ-valued. Then, for 
every t > 0, the entropy solution u(-,t) to (7) is a.e. IZ-valued. 

Proof. The proof follows from Glimm's scheme. By Lemma 2.1 and 
(ii) of Proposition 4.1, the Glimm approximation (20)-(21) is a.e. 7£-valued 
if uq is. Since 1Z is closed, Theorem 2.3 implies the result. □ 

We can now state the main result of this section: 

Proposition 5.1. Theorem 2.2 holds when the initial density profile is 
a.e. IZ-valued. 
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The proof of this proposition, from the previously established 7£-valued 
Riemann hydrodynamics (Corollary 4.2), is analogous to the proof of Theo- 
rem 3.2 on pages 16-17 of [5]. There we showed the following: for a system 
with finite initial configurations and compactly supported initial density 
profile, with the notations of Theorem 2.2 and (33), A N (rjj^ t ,u(-,t)) can- 
not grow more than o(e) on a time interval of order e. This implied the 
hydrodynamic limit for compactly supported initial profiles, and extension 
to general initial profiles followed from finite propagation arguments. The 
proof of this growth estimate for relied on macroscopic stability, and 
Lemmas 3.5 and 3.6 of [5]. Lemma 3.6 allowed us to replace u(-,t), with error 
o(e), by a piecewise constant profile with step length at least e. Lemma 3.5 
showed that, starting from this approximate profile, the hydrodynamic limit 
at times of order e was indeed the entropy solution given by Lemma 2.1 as a 
succession of noninteracting Riemann waves. In the present case, Lemma 3.5 
of Bahadoran et al. [5] is still true for 1Z- valued piecewise constant profiles: 

Lemma 5.4. Assume e > and uo(-) is an TZ-valued step function with 
step lengths at leasts; then Theorem 2.2 holds up to time min(e'/(2V),e/(2V / )), 
where V and V' are the constants from Proposition 2.1 and Lemma 5.2. 

This lemma can be derived in the same way as Lemmas 3.4 and 3.5 of 
[5], from: (a) Riemann hydrodynamics (here Corollary 4.2 instead of Theo- 
rem 2.1 in the former paper), (b) macroscopic stability, (c) finite propagation 
property for particle systems and entropy solutions (Lemma 5.2 and Propo- 
sition 2.1 above). The idea is twofold: First, Riemann hydrodynamics were 
established in Corollary 4.2 for particular random initial configurations, but 
by macroscopic stability they can be extended to any initial random se- 
quence with Riemann profile. Next, up to time min(e/(2V),e/(2V / )), we 
have a succession of noninteracting Riemann problems, both on PDE and 
particle level: this is a consequence of Lemma 5.2 and Proposition 2.1. The 
only difference here compared to [5] is that, since point (a) is only established 
for TZ- valued Riemann profiles, the conclusion also holds only for TZ- valued 
step functions uq(-). 

At this stage a new difficulty appears compared to [5]: since the above 
lemma requires TZ- valued piecewise constant profiles, we now need to ap- 
proximate the entropy solution u(-,t) — which we know is a.e. 7£-valued by 
Lemma 5.3 — with a piecewise constant profile also 7£-valued. Hence Lemma 
5.5 below is a refinement of Lemma 3.6 from [5], which requires a slightly 
different proof. 

Lemma 5.5. Assume u(-) has compact support and finite variation, and 
is a.e. TZ-valued. Let 5 > 0. Then, for e > small enough, there exists an 
approximation u £ ' S of u with the following properties: u £ ' S is a piecewise 
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constant TZ-valued function with compact support, step lengths at least e, 
and A(u e ' S , u) < e8. 

Proof. Since u has finite variation, it has left/right-hand limits u(x±0) 
at every i£K, and u(x + 0) = u(x — 0) outside an at most countable set of 
points; note that u(x ± 0) G 1Z, because 1Z is closed. Moreover, denoting by 
5(x) := \u(x + 0) — u(x — 0)| the absolute jump at x, we have J2xeR^( x ) < 
+oo. Thus there is a finite set of points J = {x±, . . . , x n }, with x\ < ■ ■ ■ < x n , 
such that 

(35) -£S(x)<6/4. 

x£J 

Let xq < x\ and x n+ \ > x n be such that the support of u is contained 
in [xo,x n+ i]. We can further divide each interval I k = (x k ,x k +\), for k = 
0, . . . ,n, into open subintervals I k) i with lengths larger than e but smaller 
than 2e. For any interval I C M, we set 

lo(I; rf) = esssup{|u(x) — u(y)\ : x, y 6 /, \x — y\ < rj}, 

uj(I) = limo;(/; n) = sup<5(x). 
vio x&i 

In particular we have sup fc=0 n uj(I k ) < 5/4. Thus, for small enough e, 

(36) (5/4 + 2 sup uj(I k ;2e) <5. 

k=0,...,n 

We shall define u £ ' S as a piecewise constant function with value outside 
(xo, x n+ i) and constant value on each interval I k i. This value can be chosen 
as follows. Set p k j = essinf/ fc ; u, p k > l = esssup /fc ; u, U (I k ,i) = {u(x ± 0), x G 
I k> l} C U. Then 

(37) Vpe(p fcjI ,p fc ' , ),3r€(«(J fc> ,)n^)Ui7(J fc ,,), \p-r\<u{I kjl )/2. 

Hence there exists r k j G 1Z such that \u k j — r k ^\ < w(I/ Cj /)/2, where u ki i de- 
notes the mean value of u on I k i. We set u £ ' 5 = r k i on I k j. 

We claim that the function u e,s thus defined satisfies the result of the 
lemma. Indeed, the integral J^—Ji^y) — u £ ' S (y)] dy has two contributions. 
The first one comes from integrating over intervals I k i that do not contain x. 
By construction, this contribution is bounded in modulus by 2eJ2k,i;xfi k t w (-^fc, 
The second contribution is the integral over part of the unique I k i contain- 
ing x, that is, f£ [u(y) — u e '^(y)]dy, where x k j is the left-hand extremity 
of Ikj. Since 

\u(x) - r k ,i\ < \u(x) -u k ,i\ + \u k ,i - r kt i\, 

the latter contribution is bounded by 2eLu{I k j] 2e) + 2euj{I k ^) /2. The result 
then follows from (35) and (36). □ 

Proceeding from Lemmas 5.3, 5.4 and 5.5 above exactly as on pages 16-17 
of [5] establishes Proposition 5.1. 
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5.2. Hydrodynamics for arbitrary entropy solutions. In order to approx- 
imate a general entropy solution with an 7?.-valued entropy solution, we use 
the following facts: 

(i) an arbitrary initial datum with compact support can be A-approximated 
with an 1Z- valued initial datum (even a {0, if}-valued initial datum is suffi- 
cient); 

(ii) the entropy solution is A-stable w.r.t. its initial datum (see [17]): 

Lemma 5.6. Assume uq(-) and vq(-) have bounded support [then, by 
Proposition 2.1, so do the corresponding entropy solutions u(-,t) and v(-,t) 
at time t]. Then A(u(-,t),v(-,t)) < A(tto(-),i>o( - ))- 

Note that this is a macroscopic analogue of Lemma 5.1. Equipped with 
the above results, we can now conclude the proof of Theorem 2.2 by deduc- 
ing general hydrodynamics from 7£-valued hydrodynamics. We shall use a 
technical result (whose proof is left to the reader): 

Lemma 5.7. Assume (r/ , N € N*) is a sequence of X-valued random 
variables such that: 

(i) a N (r] N ,dx) [defined in (IS)] converges in law along some subse- 
quence to a random measure a{dx), necessarily of the form a(dx) = u(-) dx, 
where the random density profile u(-) is a.s. a.e. [0, if]- valued. 

(ii) There exists some a € R such that, with probability tending to 1 as 
N — > oo, rj N has no particle outside the site interval [—aN,aN]. Then, for 
every v(-):R— > [0,K] with bounded support, A N (rj N ,v(-)) converges in law 
along the same subsequence to A(u(-),v(-)) . 

Proof of Theorem 2.2. In the two steps below we use the following 
notation. (j]q , iVGN*) denotes a sequence of random initial configurations 
with density profile uo(-). The random measure at(dx) =ut(-)dx (where the 
random density profile ut{-) is a.s. a.e. [0, if]-valued) denotes a subsequen- 
tial limit in law of the sequence a N (rj^ t ,dx) as N — > oo. Recall that u(-,t) 
denotes the entropy solution at time t to (7) with Cauchy datum uo(-). For 
n 6 N, Uq (•) denotes some approximation (to be defined) of uo(-) as n — > oo, 
u n (-,-) the corresponding entropy solution; (f?^'™, N € N*) a sequence (for 
fixed n ) of random initial configurations with profile Uq (■) as N — > oo, and 
r] N '™ the corresponding evolved configuration at time Nt. r/ N and tj ' n are 
coupled via the basic coupling. 

Step 1. We assume that, for some a > 0, tto(-) vanishes outside [—a, a] and 
rjQ has a.s. no particle outside [— Na, Na]. There is a sequence (uq (•), n 6 N) 
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such that Uq vanishes outside [—a, a], is a.e. {0, if }-valued, and A(«q,u ) — * 
as n — ► oo. For instance, 

n-1 
i=0 

where 5 = 2a/n, X{ = —a + id, and Oj = (KS) -1 J^ +1 uq{x) dx: indeed, since 
Uq has the same integral as uq over each interval [xi, Xj+i] , we have A(uq , uq) < 
K5 = 2aK/n. Given Uq, we can choose the sequence (r]Q ,n ,N € N) with a.s. 
no particle outside [—Na,Na\. 

Since is a.e. {0, -fTj-valued, u n (-,-) is a.e. 7£-valued by Lemma 5.3. 

Section 5.1 established the hydrodynamic limit for 7£-valued entropy solu- 
tions, hence a N (t) n '™ , dx) converges in probability to u n (-,t)dx as N — > oo 
for each n E N. Furthermore, Lemma 5.2 implies that the sequences rj^'J 1 
(for fixed n) and rj^ t satisfy assumption (ii) of Lemma 5.7. Lemma 5.7 
applied to both sequences, and the triangle inequality for A, imply that 
along some subsequence A N (rj^,rj^ t ) converges in law to A(u n (-,t),ut(-)) 
as N ^ oo. It follows from Lemma 5.1 that A(u n (-, t), u t (-)) < A(uq(-),uq(-)) 
a.s. with respect to the law of lij(-). Hence A(ii n (-, t), «t(-)) — > in proba- 
bility as n — > oo. By Lemma 5.6, A(n n (-, t), u(-, t)) — > as n — ► oo; hence 
ut{-) = u(-,t) a.s. 

Step 2. Now we make no further assumption on ito(') and • ^or n € 
we set u% = u l[-n,n] and r]^ n {x) = ^ (x)l { _ NnjNn] (x). By Lemma 5.2, r/$ t 
and r/jyj™ coincide on the site interval [iV(— n + V't),N(n — V't)] n Z with 
probability tending to 1 as iV — > oo. This and step one applied to r^J™ implies 
that ut(-) coincides with u n (-,t) a.s. on [—n + V't,n — V't]. Proposition 2.1 
implies u n (-,t) = u(-,t) a.e. on [— n + Vt, n — Vt\. Letting n — > oo yields 
■Ut(-) = u(-,t) a.s. □ 

APPENDIX A: PROOF OF PROPOSITION 2.2 

Vol'pert [26] considers the class BV of real- valued functions u(x,t) on 
R x R + * defined as follows: u £ BV if and only if, for every compact subset 
K of M x there exists a constant C(K) such that 

<<7(#)supM, 
<C(liQsupM, 

K 

for every smooth test function ip with support contained in K. Equivalently, 
u € .By if and only if its gradient in the sense of distributions is a locally 
finite vector measure. To summarize the results of Vol'pert [26] on the struc- 
ture of BV functions, we need the following notions. According to Definition 



K 



u(x, t) dt<p(x, t) dx dt 



K 



u(x, t) d x ip(x, t) dx dt 
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2.1 of [6], we say that u has an approximate jump discontinuity at y = (x,t) 
in the direction of vector v ^ 0, if there exist u~ ^ u + such that 



and v ■ z denotes scalar product in R 2 . For u = (1, — v), the meaning is that u 
has a local discontinuity at y traveling with speed v, with value u~ on the left 
and u + on the right. On the other hand, if (38) holds with u~ =u + , then it 
actually holds for any v ^ 0, and we say y is an approximate continuity point 
for u. The main results are: (i) Outside a set of one-dimensional Hausdorff 
measure 0, every point is either a point of approximate continuity or a point 
of approximate jump discontinuity, (ii) The set of approximate jump points 
is of Hausdorff dimension 1 and has locally finite one-dimensional Hausdorff 
measure. 

If u has locally bounded space variation, its spatial derivative d x u in 
distribution sense is a locally finite measure. Following (7) and Lipschitz 
continuity of G, its time derivative dtu is also a locally finite measure; hence 
u £ BV. One technical problem arises from the coexistence of two a priori 
different notions of limits for u. Indeed, (i) on the one hand we may view 
u(-,t) as a function with locally bounded space variation for all t, hence it 
has one-sided limits at every i£l for every fixed t; (ii) on the other hand, 
u is a BV function of the two variables, and as such it has approximate 
jump discontinuities as defined in (38). Limits in the sense (i) do not have 
a normal space-time vector (or equivalently a local velocity) attached to 
them, so they may not be limits in the sense (ii). Conversely, limits in the 
sense (ii) are not defined for all t and x, but only in space-time average, so 
they may not be limits in the sense (i). However (a spatially localized version 
of) Theorem 2.6 of [6] shows that these two notions coincide, provided we 
can prove the additional estimate 



for any < s < t and x,y G R, where the constant M is uniform when 
(s,t,x,y) varies in a bounded set. Equation (39) is a consequence of (7) 
and Lipschitz continuity of G. Theorem 2.6 of Bressan [6] states that: (i) 
up to a set of one-dimensional Hausdorff measure 0, the jump set of u as 
a BV function of (x,t) in the sense of Vol'pert [26] coincides with the set 
{(x,t) fix R+*; u(x — 0, t) 7^ (x + 0,t)}, (ii) the local normal is almost 
nowhere (w.r.t. one-dimensional Hausdorff measure) parallel to the time 



(38) 




where 




if v ■ z > 
if v ■ z < 



(39) 
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axis, (iii) for a jump point (x,t), let u + (resp. u~) denote the approximate 
limit on the positive (resp. negative) side of the normal chosen with positive 
x-coordinate. Then, a.e. w.r.t. one-dimensional Hausdorff measure, we have 
u ± = u(x±0,t). 

Equation (8) is satisfied by all entropy-flux pairs, if it merely holds for 
the so-called Kruzkov entropy-flux pairs ((f) c ,ip c ) defined for c £ [0,K] by 

4> c (u) := \u-c\, 4> c (u) = sgn(n - c)[G(u) - 67(c)] 

(see, e.g., Chapter 2 of [25]). The first theorem on page 256 of [26] states 
that a weak solution u £ BV of (7) satisfies (8) for (cj),ip) = (cpcipc) if and 
only if 

(40) [Mu + ) - <t>c{ur)\v t + [Mu + ) - il> B (vr)]v x < o 

holds for a.e. jump point of u w.r.t. one-dimensional Hausdorff measure, 
where u + (resp. u~) is the approximate limit of u on the positive (resp. 
negative) side of the normal vector v = {y x , vt). [Note: this is an elementary 
result if u is piecewise continuously differentiable, because then the l.h.s. of 

(40) is the trace of the l.h.s. of (8) in distribution sense along the shock line, 
when (j) = <p c and ip = ip c .) Equation (40) does not depend on the choice of 
u, but by (ii) above we may choose u with positive x-coordinate; then, by 
(iii), ti ± = u(x ±0,t). Since the l.h.s. of (40) is continuous w.r.t. c, (40) is 
satisfied simultaneously for all c outside a single exceptional set. We conclude 
by noticing that the following are equivalent for a pair (u~,u + ): 

(a) (40) holds for every c £ [0,K], 

(b) (u~,u + ) satisfies Oleinik's entropy condition (11). 

To see this, first take c = to obtain that —ft/^x is the slope of the chord 
of G between u~ and u + ; then let c vary between u~ and u + . 

APPENDIX B: PROOF OF THEOREM 2.3 

In the sequel we incorporate dependence on the discretization step Ax 
and on the sampling sequence (a) = (afc) fcgZ + £ ( — 1, l) z+ in the Glimm ap- 
proximation of u, denoting it by u Ax '( a >. We will denote by v the product 
uniform measure on ( — 1, l) z . 

For m £ N, let it mj o denote an approximation of no such that: (i) u m fi has 
locally bounded variation, (ii) u m fl converges to uq in L 1 1 oc (IR) as m—> oo. 
We denote by u m the entropy solution to (7) with initial datum u m fi, and 

by Um X '^ the Glimm approximation of u m . We know that Theorem 2.3 is 
already true when the initial datum has locally bounded variation, hence for 
u m fi. Thus, by diagonal extraction, there is a sequence 8 n | (as n — ► oo) 
and a subset A of (— 1, l) z+ with u(A) = 1, such that 

(41) V(a)£^,?n£N lim u^ x ' ia) = u m in 1^(1 X !+*). 

Ax=5„l0 
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On the other hand, by (10), u m converges to u in L\ 



*) as m 



00. 

The idea is now to compare the Glimm approximations u Ax '^ and Um X '^ ■ 
Set V" := V + R , with V the propagation speed in (10) and R the ratio 
defined in (19). We claim that, for every x < y and t < (y — x)/(2V"), 

ry-V't 



(42) E„ 



x+V"t 



\u Ax ^(z,t)-u Ax ^ a) (z,t)\dz 



< 



\u (z) -u mfi (z)\dz, 



where E^ means that (a) is integrated w.r.t. v. The meaning of (42) is that 
the Glimm approximations "nearly" satisfy the same contraction-propagation 
property (10) (with a different propagation speed V") as the real entropy 
solutions. However "nearly" means that this is not true for a fixed choice of 
the sampling sequence (a), but in average with respect to this sequence. Be- 
fore proving (42), let us conclude the proof of Theorem 2.3: for any bounded 
time interval I C (0, (y - x)/(2V")), 

rv-V"t 



E, 



x+V"t 

y-V 



u{z,t)-u Ax ^ a \z,t)\dzdt 



< 



(43) 



i" Jx+V"t 



\u(z, t) — u m (z, t) I dz dt 



y-V"t 

x+V"t 
y-V"t 



Ux+V"t 



\u m (z,t) -u Ax ' {a) {z,t)\dzdt 



u Ax '^{z,t)-u Ax '^(z,t)\dzdt 



The first term on the r.h.s. of (43) vanishes as m — > 00 by (10); the second one 
vanishes as Ax = 6 n J, for any fixed m, by (41) and dominated convergence; 
and the third one vanishes as Ax — > and then m — > 00 by (42). We have 
thus proved that, for any [x,y] C M and bounded time interval / C (0, (y — 
x)/(2V'% the limit 

r ry-V"t 

(44) lim / / \u(z,t)-u Ax ' ia) (z,t)\dzdt = 

Ax->0 J 1 J x+ V"t 

holds in L^iy) along the subsequence Ax = S n . Hence, taking growing count- 
able families of intervals [x,y] and /, and using diagonal extraction, we can 
construct a further subset A' of (— l,l) z+ with v(A') = 1, and a further 
subsequence (h n ) of (S n ), such that (44) holds along Ax = h n for every 
(a) G A', [x,y] € R and I C (0, (y - x)/{2V")). This concludes the proof of 
Theorem 2.3. 

We now proceed to prove (42). Let us introduce the following notations: 
(St)t>o denotes the evolution semigroup defined by the entropy solution to 
(7), that is, StUQ is the entropy solution at time t when the initial datum is 
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uq\ for afc £ (— 1, 1) and Ax > 0, T ak,Ax denotes the approximation operator 
denned by u+ = T ak ' Ax u k ~ in (21). Thus, for t E [t k ,t k+l ), 

(45) u Ax >( a \-,t) = S t - tk T ak ' Ax S A tT ak - lAx ■ ■■S At T ao ' Ax uo. 
It follows from definition of T ak ' Ax that 

ry+Ax 

< / \v{z) — w{z)\ dz 

J x— Ax 

for every measurable [0, K] -valued functions v(-) and w(-) on M, and every 
subinterval [x,y] Ct. On the other hand, by (10), we also have 

ry py+VAt 
(47) / \S At v(z)-S At w{z)\dz< \v(z)-w(z)\dz; 

Jx Jx-VAt 

(42) follows from (45), (46) and (47). 
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